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Numerous Earth-destroying doomsday scenarios have recently been analyzed, including break¬ 
down of a metastable vacuum state and planetary destruction triggered by a “strangelet” or micro¬ 
scopic black hole. We point out that many previous bounds on their frequency give a false sense of 
security: one cannot infer that such events are rare from the the fact that Earth has survived for 
so long, because observers are by definition in places lucky enough to have avoided destruction. We 
derive a new upper bound of one per 10 9 years (99.9% c.l.) on the exogenous terminal catastro¬ 
phe rate that is free of such selection bias, using planetary age distributions and the relatively late 


formation time of Earth. 

I. INTRODUCTION 

As if we humans did not have enough to worry about, 
scientists have recently highlighted catastrophic scenarios 
that could destroy not only our civilization, but perhaps 
even our planet or our entire observable universe. For in¬ 
stance, fears that heavy ion collisions at the Brookhaven 
Relativistic Heavy Ion Collider (RHIC) might initiate 
such a catastrophic process triggered a detailed technical 
report on the subject focusing on three risk cate¬ 
gories: 

1. Initiation of a transition to a lower vacuum state, 
which would propagate outward from its source at 
the speed of light, possibly destroying the universe 
as we know it [2,0 0 ■ 

2. Formation of a black hole or gravitational singular¬ 
ity that accretes ordinary matter, possibly destroy¬ 
ing Earth a a 

3. Formation of a stable “strangelet” that accretes or¬ 
dinary matter and converts it to strange matter, 
possibly destroying Earth id 

Other catastrophe scenarios range from uncontroversial 
to highly speculative: 

4. Massive asteroid impacts, nearby supernova explo¬ 
sions and/or gamma-ray bursts, potentially steril¬ 
izing Earth. 

5. Annihilation by a hostile space-colonizing robot 
race. 

The Brookhaven report Q concluded that if 1-3 are pos¬ 
sible, then they will with overwhelming likelihood be 
triggered not by RHIC, but by naturally occurring high- 
energy astrophysical events such as cosmic ray collisions. 
Risks 1-5 should probably all be treated as exogenous , 


i.e., uncorrelated with human activities and our techni¬ 
cal level of development. The purpose of the present pa¬ 
per is to assess the likelihood per unit time of exogenous 
catastrophic scenarios such as 1-5. 

One might think that since life here on Earth has sur¬ 
vived for nearly 4 Gyr (Gigayears), such catastrophic 
events must be extremely rare. Unfortunately, such an 
argument is flawed, giving us a false sense of security. 
It fails to take into account the observation selection ef¬ 
fect i i that precludes any observer from observing 
anything other than that their own species has survived 
up to the point where they make the observation. Even 
if the frequency of cosmic catastrophes were very high, 
we should still expect to find ourselves on a planet that 
had not yet been destroyed. The fact that we are still 
alive does not even seem to rule out the hypothesis that 
the average cosmic neighborhood is typically sterilized 
by vacuum decay, say, every 10000 years, and that our 
own planet has just been extremely lucky up until now. 
If this hypothesis were true, future prospects would be 
bleak. 

We propose a way to derive an upper bound on cosmic 
catastrophe frequency that is unbiased by such observer 
selection effects. We argue that planetary and stellar age 
distributions bound the rates of many doomsday scenar¬ 
ios, and that scenarios evading this bound (notably vac¬ 
uum decay) are instead constrained by the relatively late 
formation time of Earth. The idea is that if catastrophes 
were very frequent, then almost all intelligent civiliza¬ 
tions would arise much earlier than ours did. Using data 
on planet formation rates, it is possible to calculate the 
distribution of birth dates for intelligent species under 
different assumptions about the rate of cosmic steriliza¬ 
tion. Combining this with the information about our own 
temporal location enables us to conclude that the cosmic 
sterilization rate for a habitable planet is at most of order 
one per Gigayear. 
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FIG. 1: The left panel shows the probability distribution for observed planet formation time assuming catastrophe timescales r of oo 
(shaded), 10, 9, 8, 7, 6, 5, 4, 3, 2 and 1 Gyr, respectively (from right to left). The right panel shows the probability of observing a 
formation time >9.1 Gyr (that for Earth), i.e., the area to the right of the dotted line in the left panel. 


II. AN UPPER BOUND ON THE 
CATASTROPHE RATE 

Suppose planets get randomly sterilized or destroyed at 
some rate r _1 which we will now constrain. This means 
that the probability of a planet surviving a time t decays 
exponentially, as e -t / r . 

The most straightforward way of eliminating observer 
selection bias is to use only information from objects 
whose destruction would not yet have affected life on 
Earth. We know that no planets from Mercury to Nep¬ 
tune in our solar system have been converted to black 
holes or blobs of strange matter during the past 4.6 
Gyr, since their masses would still be detectable via their 
gravitational perturbations of the orbits of other planets. 
This implies that the destruction timescale r must be cor¬ 
respondingly large — unless their destruction is be linked 
to ours, either by a common cause or by their implosion 
resulting in the emission of doomsday particles like black 
holes or strangelets that would in turn destroy Earth. 
This observer selection effect loophole is tightened if we 
consider extrasolar planets that have been seen to par¬ 
tially eclipse their parent star Q and are therefore known 
not to have imploded. The doomsday particles discussed 
in the literature would be more readily captured gravi¬ 
tationally by a star than by a planet, in which case the 
observed abundance of very old (> 10 Gyr) stars (e.g., 
0) would further sharpen the lower bound on r. 

The one disaster scenario that exploits the remaining 
observer bias loophole and evades all these constraints is 
vacuum decay, either spontaneous or triggered by a high- 
energy event. Since the bubble of destruction expands 
with the speed of light, we are prevented from observing 
the destruction of other objects: we only see their de¬ 
struction at the instant when we ourselves get destroyed. 


In contrast, if scenarios 2 or 3 involved doomsday particle 
emission and proceed as a chain reaction spreading sub- 
luminally, we would observe spherical dark regions cre¬ 
ated by expanding destruction fronts that have not yet 
reached us. We will now show that the vacuum decay 
timescale can be bounded by a different argument. 

The formation rate f P (t p ) of habitable planets as a 
function of time since the Big Bang is shown in Figure ^ 
(lef t panel, shaded distribution). This estimate is from 
Hoi , based on simulations including the effects of heavy 
element buildup, supernova explosions and gamma-ray 
bursts. If regions of space get randomly sterilized or de¬ 
stroyed at a rate t , then the probability that a ran¬ 
dom spatial region remains unscathed decays as e -t / T . 
This implies that the conditional probability distribu¬ 
tion f*{tp) for the planet formation time t p seen by an 
observer is simply the shaded distribution f p (t p ) multi¬ 
plied by e _tp / r an d rescaled to integrate to unity, giving 
the additional curves in Figure Q (left panel). 1 As we 
lower the catastrophe timescale r, the resulting distribu¬ 
tions (left panel) are seen to peak further to the left and 


1 Proof: Let / 0 (t 0 ) denote the probability distribution for the time 
t a after planet formation when an observer measures t p . In our 
case, t 0 = 4.6 Gyr. We obviously know very little about this 
function f 0 , but it fortunately drops out of our calculation. The 
conditional probability distribution for t p , marginalized over f 0 , 
is 

/* ^ tp 

f*(t P )oz fo(to)fp(t p )e — dt 0 oc fp(t p )e~~, (1) 
J 0 

independently of the unknown distribution f 0 (t 0 ), since 
e~(t 0 +tp)/T — e —t 0 lT e —tplT anc j ^ 0^0 the entire integrand is 
separable into a factor depending on t p and a factor depending 
on t Q . 
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the probability that Earth formed as late as observed (9.1 
Gyr after the Big Bang) or later drops (right panel). The 
dotted lines show that we can rule out the hypothesis that 
r < 2.5 Gyr at 95% confidence, and that the correspond¬ 
ing 99% and 99.9% confidence limits are r > 1.6 Gyr and 
r > 1.1 Gyr, respectively. 

Risk category 4 is unique in that we have good direct 
measurements of the frequency of impacts, supernovae 
and gamma-ray bursts that are free from observer se¬ 
lection effects Our analysis therefore used the habitable 
planet statistics from [1C| that folded in such category 4 
measurements. 

Our bounds do not apply in general to disasters of 
anthropogenic origin, such as ones that become possible 
only after certain technologies have been developed, e.g., 
nuclear annihilation or extinction via engineered microor¬ 
ganisms or nanotechnology. Nor do they apply to natu¬ 
ral catastrophes that would not permanently destroy or 
sterilize a planet. In other words, we still have plenty to 
worry about pmmm However, our bounds do ap¬ 
ply to exogenous catastrophes (e.g., spontaneous or cos¬ 
mic ray triggered ones) whose frequency is uncorrelated 
with human activities, as long as they cause permanent 
sterilization. 

Our numerical calculations made a number of assump¬ 
tions. For instance, we treated the exogenous catastrophe 
rate r _1 as constant, even though one could easily imag¬ 
ine it varying by of order 10% over the relevant timescale, 
since our bound on r is about 10% of the age of the Uni¬ 
verse. 2 Second, the habitable planet formation rate in¬ 
volved several assumptions detailed in jl(j which could 
readily modulate the results by 20%. Third, the risk 
from events triggered by cosmic rays will vary slightly 
with location if the cosmic ray rate does. Fourth, due 
to cosmological mass density fluctuations, the mass to 
scatter off of varies by about 10% from one region of 
size cr ~ 10 9 lightyear region to another, so the risk of 
cosmic-ray triggered vacuum decay will vary on the same 
order. 

In summary, although a more detailed calculation 
could change the quantitative bounds by a factor of order 
unity, our basic result that the exogenous extinction rate 
is tiny on human and even geological timescales appears 


rather robust. 


III. CONCLUSIONS 

We have shown that life on our planet is highly unlikely 
to be annihilated by an exogenous catastrophe during 
the next 10 9 years. This numerical limit comes from the 
scenario on on which we have the weakest constraints: 
vacuum decay, constrained only by the relatively late for¬ 
mation time of Earth, conclusion also translates into a 
bound on hypothetical anthropogenic disasters caused by 
high-energy particle accelerators (risks 1-3). 

This holds because the occurrence of exogenous catas¬ 
trophes, e.g., resulting from cosmic ray collisions, places 
an upper bound on the frequency of their anthropogenic 
counterparts. Hence our result closes the logical loop¬ 
hole of selection bias and gives reassurance that the risk 
of accelerator-triggered doomsday is extremely small, so 
long as events equivalent to those in our experiments oc¬ 
cur more frequently in the natural environment. Specif¬ 
ically, the Brookhaven Report Q suggests that possible 
disasters would be triggered at a rate that is at the very 
least 10 3 times higher for naturally occurring events than 
for high-energy particle accelerators. Assuming that this 
is correct, our 1 Gyr limit therefore translates into a con¬ 
servative upper bound of 1/10 3 x 10 9 = 10 -12 on the an¬ 
nual risk from accelerators, which is reassuringly small. 


2 As pointed out by Jordi Miralda-Escude (private communica¬ 
tion) , the constraint from vacuum decay triggered by bubble nu- 
cleation is even stronger than our conservative estimate. The 
probability that a given point is not in a decayed domain at time 
t is the probability of no bubble nucleations in its backward light 
cone, whose spacetime 4-volume tx //' for both matter-dominated 
and radiation-dominated expansion. A constant nucleation rate 
per unit volume per unit time therefore gives a survival prob¬ 
ability e~ b/ T ) for some destruction timescale r. Repeating 
our analysis with e -t / T replaced by the sharper cutoff e — d/' r ) 
sharpens our constraint. Our quoted bound corresponds to the 
conservative case where r greatly exceeds the age of the universe 
at the dark energy domination epoch, which gives a backward 
lightcone volume oc t. 
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